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Рассмотрен новый подход к решению задачи переноса излучения в дисперсной среде, базирующийся на использовании инди- 
катрисы рассеяния излучения в виде интегральных параметров, вычисленных на основе двух индикатрис, представленных во 
взаимно перпендикулярных плоскостях. Такое представление индикатрисы рассеяния излучения адаптировано к потоковым ме- 
тодам решения задачи переноса излучения. Анализ результатов решения производится с помощью коэффициентов асимметрии 
индикатрисы рассеяния излучения и тела яркости рассеивающего объема. 


Вопрос о прохождении излучения через диспе- 
рсную среду, состоящую из частиц неправильной 
формы, является одним из актуальных в теории пе- 
реноса излучения. Представление индикатрисы 
рассеяния излучения в этом случае аналогично 
представлению индикатрисы рассеяния для частиц 
правильной формы, а для системы частиц непра- 
вильной формы с произвольной ориентацией в 
пространстве используется усредненная индикат- 
риса рассеяния излучения. 

Стандартным условием в теории рассеяния из- 
лучения является условие осесимметричности ин- 
дикатрисы рассеяния излучения. Именно при этом 
условии определяются некоторые интегральные 
параметры индикатрисы рассеяния, например, ко- 
эффициент асимметрии индикатрисы рассеяния 
излучения или средний косинус угла рассеяния [1]. 
Для характеристики степени анизотропии рассея- 
ния часто этих параметров бывает достаточно, осо- 
бенно когда требуется оперативность в получении 
информации об угловом распределении рассеян- 
ного излучения, хотя иногда возникает неодноз- 
начность, например, коэффициенты асимметрии 
релеевской индикатрисы и сферической одинако- 
вы. В случае рассеяния излучения на частице неп- 
равильной формы для характеристики анизотро- 
пии рассеяния одного коэффициента асимметрии 
индикатрисы рассеяния, характеризующего рассе- 
яние в переднюю и заднюю полусферы, недоста- 
точно. Требует уточнения в этом случае и стандарт- 
ное представление индикатрисы рассеяния излуче- 
ния, и методики ее получения. 

Для частиц неправильной формы положение о 
симметрии рассеяния не выполняется, и, соответ- 
ственно, интегральные характеристики рассеяния, 
введенные для симметричной индикатрисы рассе- 
яния излучения, не отражают особенностей рассе- 
яния излучения на таких частицах. Более того, 
обычная "плоская" индикатриса рассеяния излуче- 
ния не дает полной информации о пространствен- 
ном угловом распределении излучения при рассея- 
нии на частицах неправильной формы. В общем 
случае полную информацию об угловом распреде- 
лении излучения, рассеянного частицей непра- 
вильной формы, может дать только объемная ин- 
дикатриса рассеяния излучения. Получение такой 
индикатрисы рассеяния теоретически, а особенно 


экспериментально, весьма трудоемко. Кроме того, 
индикатриса рассеяния излучения является одной 
из основных характеристик, определяющих пере- 
нос излучения в дисперсной среде, поэтому форма 
ее представления должна быть адаптирована к ме- 
тоду, используемому для расчета характеристик пе- 
реноса излучения [2, 3]. 

В работе предлагается следующая методика 
представления индикатрисы рассеяния излучения 
при рассеянии на частицах неправильной формы. 
Рассмотрим рассеяние излучения на частице неп- 
равильной формы в декартовой системе координат, 
когда излучение падает в направлении оси х. Объ- 
емная индикатриса рассеяния излучения заменяет- 
ся двумя индикатрисами рассеяния, определенны- 
ми в двух взаимно перпендикулярных плоскостях 
ху и хі. Обозначим их Хх; ( Ѳ) и %„{ Ѳ) соответственно, 
где Ѳ - угол рассеяния. Интегральные параметры 
определяются в общем случае следующим образом 

Иу = | Х ху (Ѳ)сЮ.-, 

(П) 

Я- = I Х Х ;(Ѳ )сЮ., 

(П) 

где О - телесный угол. 

Подробно все интегральные параметры по осям 
декартовой системы координат для двух индикат- 
рис рассеяния излучения определяются следую- 
щим образом 
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Определим коэффициенты асимметрии по 
осям координат для двух индикатрис рассеяния 

Рі + /У (Р* + Р\ + Рз + Р У 6 ) 
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Предложенные коэффициенты позволяют 
учесть асимметрию рассеяния по осям координат 
и, тем самым, устранить неоднозначность связи 


между параметрами среды и излучения. Данное 
обстоятельство особенно важно в случае частиц 
неправильной формы, в силу различных физичес- 
ких причин ориентированных в пространстве, нап- 
ример, в конвективных потоках, электрических по- 
лях или в сплошных средах. 

Индикатриса рассеяния излучения определяет 
угловое распределение рассеянного излучения для 
одной частицы, или микрообъема, в котором про- 
исходит только однократное рассеяние. Угловое 
распределение излучения с учетом многократного 
рассеяния внутри дисперсной среды характеризу- 
ется полем или телом яркости. Однако, определе- 
ние поля яркости внутри среды в ряде случаев свя- 
зано со значительными трудностями, поэтому 
обычно рассматривают угловое распределение из- 
лучения вне рассеивающего объема. По аналогии с 
коэффициентами асимметрии индикатрисы рассе- 
яния целесообразно ввести коэффициенты асим- 
метрии излучения, рассеянного макрообъемом. 
Аналогично формулам для определения коэффи- 
циентов асимметрии индикатрисы рассеяния, ани- 
зотропию рассеяния излучения макрообъема будем 
характеризовать следующими коэффициентами 
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Здесь І +х - поток излучения, выходящий из сре- 
ды в направлении падающего потока вдоль оси х, 
/ х - поток излучения, отраженный дисперсной 
средой, Р у , Р, - потоки рассеянного излучения, вы- 
ходящие из среды в направлениях +у, +х 

Примером использования вышеприведенных 
формул являются результаты расчетов (рис. 1, 2), 
позволяющие проследить деформацию тела яркости. 


Аз 



дольной оптической толщины дисперсной среды: 
1, 3) Ту = 1; 2, 4) Ту = 100 ; 1,2) а=1 ; 3, 4) О 


Рассмотрим результаты исследований формы 
тела яркости пространственно ограниченного рас- 
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сеивающего объема в зависимости от оптических 
размеров и параметров среды. Форма тела яркости 
объема оценивалась с помощью формул (*). 

На рис. 1 представлены результаты расчетов ве- 
личины коэффициента асимметрии Аз* в зависи- 
мости от оптической толщины дисперсной среды 
для сферической (а=1) и анизотропной (С1) [4] ин- 
дикатрис рассеяния излучения. 

Расчеты проводились по методике, изложенной 
в [5], для двух индикатрис рассеяния излучения, ха- 
рактеризуемых коэффициентами асимметрии а й = 1 
(сферическая индикатриса) и а л ~\2 (индикатриса 
С1), при поперечных оптических размерах т у =т=1, 
10, 100 и вероятности выживания кванта Л=1. Оче- 
видно, что при таких исходных данных тело яркос- 
ти имеет осесимметричное строение. По мере уве- 
личения оптических размеров среды доля многок- 
ратно рассеянного излучения в общем потоке воз- 
растает, в то время как доля излучения, прошедше- 
го среду без рассеяния и ослабленного по экспо- 
ненциальному закону, уменьшается. Этим фактом 
объясняется более резкая зависимость коэффици- 
ента А? от оптической толщи г для среды, характе- 
ризуемой анизотропной индикатрисой рассеяния 
излучения. Следует отметить, что при всех оптичес- 
ких поперечниках, использованных в расчетах, чет- 
ко определяется граничное значение продольной 
оптической плотности дисперсной среды, выше ко- 
торой коэффициент асимметрии не меняется, а, 
следовательно, не меняется и тело яркости рассеи- 
вающего объема. Неизменность тела яркости при 
увеличении оптической глубины рассеивающей 
среды обычно характеризует наступление глубин- 
ного режима. В данном случае удивляет достаточно 
раннее наступление глубинного режима (ту-10), что 
может объясняться условиями освещения дисперс- 
ной среды (бесконечно широкий пучок) [3]. 

Расчеты зависимости величины коэффициента 
асимметрии от поперечных оптических размеров 
дисперсной среды выполнены для тех же двух ин- 
дикатрис рассеяния излучения, что и в предыду- 
щем случае, и консервативной среды. Полученные 
данные показывают достаточно слабую зависи- 
мость коэффициента асимметрии А ? от попереч- 
ных оптических размеров дисперсной среды. Это 
происходит даже в том случае, если изменяются 
оба поперечных оптических размера среды. В то же 
время сохраняется сильная зависимость коэффи- 
циента асимметрии от вида индикатрисы рассея- 
ния излучения. Исследование зависимости Аз от 
поперечных оптических размеров позволяет опре- 
делить предельные оптические размеры, при кото- 
рых среду можно считать пространственно неогра- 
ниченной в поперечном оптическом сечении. Пос- 
тоянство коэффициента асимметрии при измене- 
нии оптических размеров дисперсной среды гово- 
рит также о сформировавшемся теле яркости, что 
аналогично наступлению глубинного режима. 

Зависимость деформации тела ярости от наличия 
поглощения в дисперсной среде приведена на рис. 2. 


Кривые 1, 2 получены для сферической индикатри- 
сы рассеяния, кривые 3, 4 - для индикатрисы С1 [4]. 
Как видно из рисунка, влияние поглощения на фор- 
му тела яркости имеет неоднозначный характер. В 
случае дисперсной среды, имеющей форму куба с 
оптическими размерами т =2, тело яркости с увели- 
чением поглощения становится более вытянутым 
вне зависимости от формы индикатрисы рассеяния 
излучения (кривые 1, 3). С увеличением оптических 
размеров дисперсной среды зависимость формы тела 
яркости от увеличения поглощения в среде меняется 
на противоположную (кривые 2, 4). Это означает, что 
имеется дисперсная среда таких оптических разме- 
ров, для которой изменение поглощения в среде 
почти не влияет на форму тела яркости. 


Аз 45 



Рис. 2. Зависимость коэффициента асимметрии Аз от веро- 
ятности выживания кванта А: 1, 3) т = 2; 2, 4) т = 20; 
1, 2) а=1; 3, 4) О 

Таким образом, в данной работе: 

- предложен новый подход к решению задачи пе- 
реноса излучения в дисперсной среде с частица- 
ми неправильной формы, с использованием 
индикатрисы рассеяния излучения в виде ин- 
тегральных параметров, вычисленных на осно- 
ве двух индикатрис, представленных во взаим- 
но перпендикулярных плоскостях; 

- представление индикатрисы рассеяния излуче- 
ния в виде интегральных параметров адаптиро- 
вано к потоковым методам решения задачи пе- 
реноса излучения; 

- анализ результатов решения производится с по- 
мощью коэффициентов асимметрии индикат- 
рисы рассеяния излучения и тела яркости рас- 
сеивающего объема; 

- показано, что предложенные коэффициенты 
асимметрии являются информативными и 
чувствительными характеристиками, позволя- 
ющими оценить пространственное распределе- 
ние рассеянной объемом радиации. 

Изучение поведения данных коэффициентов в 

зависимости от параметров дисперсной среды поз- 
воляет оценить границы применимости экспонен- 
циального закона ослабления излучения, наступ- 
ление глубинного режима и предельные оптичес- 
кие размеры, при которых среду можно считать не- 
ограниченной. 
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